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We study the cosmology of a generalized Galileon field <j> with five covariant Lagrangians in which 
(j> is replaced by general scalar functions fi((f>) (i = 1, • ■ • ,5). For these theories, the equations of 
motion remain at second-order in time derivatives. We restrict the functional forms of fi(4>) from 
the demand to obtain de Sitter solutions responsible for dark energy. There are two possible choices 
for power- law functions fi((j>), depending on whether the coupling F(<j}) with the Ricci scalar R is 
independent of 4> or depends on cp. The former corresponds to the covariant Galileon theory that 
respects the Galilean symmetry in the Minkowski space-time. For generalized Galileon theories we 
derive the conditions for the avoidance of ghosts and Laplacian instabilities associated with scalar 
and tensor perturbations as well as the condition for the stability of de Sitter solutions. We also 
carry out detailed analytic and numerical study for the cosmological dynamics in those theories. 
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t—i • I. INTRODUCTION 

o : 

, The A-Cold-Dark-Matter (ACDM) model has been consistent with observational data, but the energy scale of dark 
, > ■ energy is too low to be compatible with the cosmological constant originated from the vacuum energy in quantum 
field theory [1]. Since the observations allow the variation of the dark energy equation of state [2], many models have 
been proposed to explain the present accelerated expansion of the Universe [3]. For example, a light scalar field with 
a slowly varying potential, called quintessence, was introduced as an alternative to the cosmological constant [4]. In 
general, however, it is not easy to construct viable particle physics models of quintessence because of an extremely 
light mass required for the cosmic acceleration today [5]. 

Another approach for addressing the dark energy problem is to modify the law of gravity from General Relativity at 

■ large distances [6]. In this approach there have been two main streams. The first consists of introducing a Lagrangian 
■ \ for gravity built up out of the Ricci. Riemann, and metric tensors, which generally leads to 4-th order differential 

^j*- equations. The f(R) gravity [7] and the Gauss-Bonnet gravity [8] belong to this class. The second consists of higher 
(-h ' dimensional models that realize the cosmic acceleration through the gravitational leakage to extra dimensions. The 
i— i Dvali-Gabadadze-Porrati (DGP) braneworld model [9] belongs to this class (see Refs. [10]). 

In general, modified gravity models of dark energy need to be constructed to recover the General Relativistic 
, behavior in the regions of high density for the consistency with local gravity experiments. In f(R) gravity, there have 

■ been a number of viable models in which a scalar-field degree of freedom ("scalaron" [11]) has a large mass in the 
] region where the Ricci scalar R is much larger than its cosmological value i?o today [12]. Provided that the chameleon 

£NJ ■ mechanism [13] is at work in the local regime, the gravitational coupling with non-relativistic matter can be suppressed 
to be compatible with solar system experiments [14]. There is also another mechanism called the Vainshtein screen 

■ effect [15] in which non-linear effects can effectively decouple the scalar field from gravity. Originally the Vainshtein 
mechanism was applied to the theories of massive gravity like Fierz-Pauli gravity [16] (see also Refs. [17]), but the 

CD non-linearities imply the presence of a ghost state in such theories [18]. 

In the DGP model non-linear field self-interacting Lagrangians such as □0(9 M c/><9 At 0) arise from a brane-bending 
mode (i.e. a longitudinal graviton) [19]. This allows the decoupling of (j> from gravitational dynamics in the local 
region. Unfortunately the self-accelerating solution in the DGP model contains a ghost mode [20] even in the absence 
of non-linear terms. Moreover the model is disfavored from the combined data analysis of supernovae la and baryon 
acoustic oscillations [21]. 

Mostly inspired by the DGP model, Nicolis et al. [22] derived the five Lagrangians that lead to the field equations 
invariant under the Galilean symmetry d^cf) —> d^+b^ in the Minkowski space-time [including the term □ 0(9^ <pd^(f))]. 
The scalar field that respects the Galilean symmetry is dubbed "Galileon". Each of the five terms only leads to second- 
order differential equations, keeping the theory free from unstable spin-2 ghost degrees of freedom. If we extend the 
analysis in Ref. [22] to the curved space-time, the Lagrangians need to be promoted to the covariant forms. Deffayet 
et al. [23, 24] derived the covariant Lagrangians d (i = 1, ■ • ■ ,5) that keep the field equations up to second-order. We 
caution that in the curved space-time the Galilean symmetry is in general broken for non-linear field self-interacting 
terms, but in the Minkowski space-time it is preserved for the covariant Lagrangians Ci (i = 1, • • • ,5). 

The (modified) Galileon gravity has been extensively applied to cosmology recently [25-41]. One application is to 
introduce the non-linear field self-interaction of the form £((/))D(/)(cL,0d A1 </)) in the action of (generalized) Brans-Dicke 
theories [26-28, 31, 32], where £ is a function of (j>. Although such a term breaks the Galilean symmetry, the field 
equations remain at second-order. Moreover, for suitable choices of the function £(0), there exist de Sitter (dS) 
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solutions responsible for dark energy even in the absence of the field potential. The presence of the non-linear term 
also allows the decoupling of the field from gravity in the regions of high density under the Vainshtein mechanism. 

Another application of Galileon gravity to cosmology is to study the expansion history of the Universe in the 
presence of the covariant Lagrangians £^ (i = 1, • • • ,5) mentioned above. The cosmology up to the term £4 has been 
discussed in Ref. [30], which showed the existence of stable dS solutions. Recently the full cosmological dynamics 
including the term £5 have been studied in Ref. [34]. The viable model parameter space has been found by studying 
the conditions for the avoidance of ghosts and Laplacian instabilities. Interestingly there exists a tracker solution 
that finally approaches a stable dS solution. The equation of state of dark energy exhibits a peculiar phantom-like 
behavior along the tracker. 

In this paper we shall study general Galileon theories in which the field 4> in the covariant Lagrangians Ci (i = 
1, ■ ■ • ,5) is replaced by general functions fi(4>). Since /i(</>) are scalar functions, the resulting field equations also 
remain at second-order. We constrain the forms of fi{4>) from the requirement to obtain dS solutions. This constraint 
gives rise to the Galileon theory with fi(<fi) cx (j> as a specific case. For general functions fi((f>) we a l so derive the 
conditions for the avoidance of ghosts and Laplacian instabilities. This is useful to constrain the viable parameter 
space of those theories. We shall perform detailed analytic and numerical study for the cosmological dynamics of 
generalized Galileon theory with several different choices of /i (</>). 



II. GENERALIZED GALILEON THEORIES 



In the curved space-time the Galilean symmetry is broken even for the Lagrangian £2 = (V</>) 2 = d^tfid^^. Then 
this symmetry is not restrictive when we study the covariant generalization of the Galileon field. On the other hand, 
the covariant Galileon formalism leads to second-order field equations. We study a Lagrangian that gives second-order 
equations of motion, such that the theories recover the covariant Lagrangian in Refs. [23, 24] as a specific case. We 
will consider two generalizations of the covariant Galileon theory: (i) scalar couplings with both the Ricci scalar R 
and the Gauss-Bonnet (GB) term Q are introduced, (ii) the covariant Galileon terms are extended to more general 
functions. 

As for the first point, this step is compatible with the approach of field theory, because such scalar couplings 
generally exist and they are consistent with general covariance (and even with the Galileon symmetry, as in the 
Minkowski background their contributions to the equations of motion of the field identically vanish) . Moreover the 
scalar couplings with R and Q give only second-order contributions. It is true that the GB term can change the Ultra- 
violet behavior for the modes, but this property also holds for all the remaining terms coming from the extended 
Galileon action. 

As for the second point, we can replace the scalar field <fi in each Lagrangian term with a function of the field 
itself. The Lagrangian £2, for example, can be modified to (V/2(</>)) 2 = d^f^cj)) 8^/2 (</>)■ The equations still remain 
at second-order because f2 is a scalar quantity itself. We will consider this generalization for all the Galileon terms, 
introducing different functions fi(4>) (i = 1,2,*--) for each of them. 

According to the above prescription, we introduce the following Lagrangians as the generalization of those introduced 
by Deffayet et al. [23]: 

Ci = h{4>), (i) 

£2 = (V/ 2 (0)) 2 , (2) 
£ 3 = (n/3(0))(V/ 3 (0)) 2 , (3) 
£4 = (V/ 4 (0)) 2 [2(D/ 4 (^)) 2 - 2/4(^/4 - (R/2) (V/ 4 (0)) 2 ] , (4) 

£5 = (V/ 5 (0)) 2 [(Q/sW) 3 - mhmhW^hW'^ +2f 5 (4>); ti U h(4>);u P m);p fl 

-^h{4>)-M^ u hW' p G vp ], (5) 

where R is the Ricci scalar and G vp is the Einstein tensor. One can also introduce the following terms 

£ 6 = F(4)R, (6) 
£7 = ^)5, (7) 

which vanish in the Minkowski space-time. Here Q = R 2 — AR a pR a P + R a p 1 sR al3 ~ >S represents the GB combination, 
where R a p is the Ricci tensor and R a p~fS is the Riemann tensor. 
The covariant action we shall discuss is given by 

S=\ f d A x^C+ f A A xC M , (8) 
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where g is a determinant of the space-time metric g^, Cm is a matter Lagrangian, and 

C = F(0)P + i{4>)G + Ci+ e 2 C 2 + C 3 + e 4 C 4 + C 5 ■ (9) 

In order to control the signs of £2,4 we have introduced the factors 62,4, which are either +1 or — 1. For the terms 
£1,3,5 we g e t the opposite signs by replacing /i(0) with —fi{<j>). For the matter Lagrangian Cm we consider the 
contribution of two perfect fluids C^l (i = 1,2), described by the barotropic equations of state of the form Wi = Pi/ pi 
(i = 1,2). Note that Pj and pi are the pressure and the energy density of fluids, with the energy-momentum tensor 

T# = -(2/V=s)AC$/*<r. 



III. BACKGROUND COSMOLOGY 

Consider the flat Friedmann-Lemaitre-Robertson- Walker (FLRW) space-time with the line element 

ds 2 = g^dx^dx" = -dt 2 + a 2 (t) dx 2 , (10) 

where a(t) is the scale factor with cosmic time t. Varying the action (8) with respect to we obtain the following 
equations of motion 

3FH 2 + 3HF + A/2 + e 2 / 2 2 /2 - 3-ff/f + Ahe 4 H 2 ff/2 - 21P 3 /| + 12# 3 £ = Pi + P2 , (11) 
3FP 2 + F + 2HF + 2FH + A/2 - e 2 f 2 /2 - f 2 f 3 + 3e 4 / 4 3 (8P/ 4 + 3P 2 / 4 + 2P/ 4 )/2 
-3P/ 5 4 (5P/ 5 + 2P/ 5 + 2P 2 / 5 ) + 4H[Hii+2(H 2 + H)£] = -{w lPl + w 2 p 2 ) , (12) 

where a dot represents a derivative with respect to t. The matter fluids obey the continuity equations 

Pi + 3H(1 + Wi)pi = 0, (* = 1,2). (13) 

Differentiating Eq. (11) in terms of t and eliminating the terms w\pi + w 2 p 2 from Eq. (12), we get the generalized 
Klein-Gordon equation for the scalar field. For the perfect fluids we consider radiation (wi = 1/3, p\ = p r ) and 
non-relativistic matter (102 = 0, p 2 = p m )- 

Let us restrict the functional forms of fi (</)), F{<f), and by demanding the existence of dS solutions responsible 
for dark energy. We shall focus on the theories in which the late-time cosmic acceleration can be realized by the field 
kinetic terms rather than the field potential, so that we set 

/i(0)=O. (14) 

The absence of this term and more in general of a potential for the field can be implemented by invoking an additional 
shift symmetry — > + c for the Galileon Lagrangian. The condition (14) is also important for another reason. 
In the Minkowski space-time (H = 0) the only solution to the equations of motion without matter corresponds to 
= 0, provided that f 2 ^ ^ 0. This implies that the field is frozen, and thus implementing the Vainshtein mechanism. 
Moreover the field perturbations would propagate with the speed of light in the Minkowski background. 
The Friedmann equation (11) can be written in the form 

ci r + n m + n U E = 1 , (15) 

where n r = p r /(3FH 2 ), VL m = p m /(3FH 2 ), and 

f fi fi 15/I 7if/f mi 

n ^ = ~HF~ €2 6FIP + HF' €i ^F~ + ~ F~ ' (16) 

In order to realize the late-time dS solutions we take the power-law functions for F{<f>), /i(0), and £(0) in terms of 0. 
We can classify the theories into two classes: (i) F is constant, and (ii) F depends on the field 0. 



A. Constant F 



At the dS point (H = Pds = constant) we require that each term in Eq. (16) does not vary in time. For constant 
F, the functions /, and £ need to have the field dependence proportional to 0. Then the dS solutions can be realized 



4 



for constant <fi. In order to have dimensionless couplings di of the orders of unity, we write the functions F, /,-, and £ 
in the forms 

where M p \ = 2.43 x 10 18 GeV is the reduced Planck mass. Together with the dimensionless constants di (i = 2, • • • , 5) 
and we have introduced the mass scales 

M = ^y /3 ff dS «10- 40 M pl , (18) 
M s = ^«10- 120 M pl , (19) 

where we have used i?ds ~ 10 _60 M p i. 

Defining the mass scales M and as given above, we can express the density parameter J7de in the form 

^de = + dixy - ^-e 4 d 4 4 y 2 +7d 5 5 ^-- id^ , (20) 

D Zi XX 

where 

<t> 2 H 2 

x = , y = x — s- . (21) 

We shall consider the case in which the dimensional variables x and y are not much different from the orders of unity 
today. In the asymptotic past we require that x -C 1, y *C 1, and y 3> x 2 to recover the General Relativistic behavior. 
If the coefficients di and d^ are of the orders of unity, the dominant contribution in Eq. (20) comes from the last term. 
In such a case, however, the GB term has been dominant during the whole cosmological evolution by today. In order 
to avoid this behavior, we set 

de = o, (22) 

when we discuss the cosmological dynamics. 

Using Eqs. (11) and (12), we obtain the following relations: 

e 2 (ta3) 2 = 6 + 9a -12/?, (23) 
(d 3 x ds ) 3 = 2 + 9a -9/3, (24) 

where x d g is the value of x at the dS point, and 

a = e 4 (cLiXds) 4 , P = (d 5 x d s) 5 ■ (25) 

We note that the theory with the functions (17) corresponds to the covariant Galileon theory discussed in Refs. [23, 31]. 



B. Non-constant F 



Let us consider theories in which F depends on the field </>. If we take the power-law function of the form F oc <fi p 
with constant p, it follows from Eq. (16) that the function <fr/(H(f>) is required to be constant at the dS point. From 
this demand we can restrict the functions F, /j, £ in the forms 

f = mj-v, h = d 2 M^- p/2 r /2 , h = d 3 A/ 3 - p/ V /3 , u = d,M- i/2 - p/ u p, \ 

h = d 5 M- i/5 - p/5 ^ , £ - d s M~ P r , (26) 
where di (i = 2, • • • , 5), d^ are dimensionless constants, Mi (i = 2, • • • ,5) and are mass scales defined by 

M \ (P-2J/P / M N (p-2)/( P +2) , M \ ( P -2)/(p+4) 

P 1 I tr._ if . — I P 1 I ej\_ ir - I P 1 



M p i, M 3 ee ( j^j H dS , M 4 = \-^\ H dS , M 5 = (-^ I H dS , M £ = M 3 . 

(27) 
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We note that there are other possibilities to obtain dS solutions, depending on the forms of F(cf>). If we choose the 
function F(<j>) oc e^, where fj, is a constant, then it is possible to realize the dS solution for the choice fa oc e^/ 2 , 
/3 oc e M< ^/ 3 , /4 oc e M< ^ 4 , and /s oc e^/ 5 . In the sense that </> = constant along the dS solution this theory is related 
with the constant F theory given in (17). Since we want to consider the case in which the dS solution is realized in a 
different way, we shall study the cosmological dynamics for the theories with (26) in which cb/4> = constant at the de 
Sitter solution. 

For the theories (26) the Galileon symmetry is explicitly broken. However, even for the Galileon action (17), the 
Galileon symmetry is restored only on the Minskowski background. Therefore the Galileon symmetry does not restrict 
the form of the Lagrangian on curved backgrounds. Nonetheless, these theories may represent an effective action for 
some more fundamental theory, e.g., extra-dimensional models. In fact the covariant Galileon is the generalization 
of the decoupling limit of the DGP braneworld model. The existence of dS solutions in this theory opens up the 
possibility of studying these generalizations of the original Galileon field. Moreover the model is not plagued by the 
Ostrogradski instability because the field equations remain at second order. The situation here is not very different 
from other dark energy models, such as f(R) gravity, where the Lagrangian is constructed by hands to realize the late- 
time cosmic acceleration but it is supposed to be originated from some fundamental theory. In this directions other 
papers appeared which tried to generalize the Galileon Lagrangian without imposing the original Galileon symmetry 
[26-28, 31, 37], as it is violated, by construction, on curved backgrounds. 

The density parameter f^DE can be expressed as 

tt DE = -Pi - ^e 2 d 2 2 p 2 x 2 + ^d 3 3P 3 xf - ^4p 4 f + ^4p 5 ^ - , (28) 

where 

6 H 

XS W '"5=- <29) 

For p and dj (i = 2, ■ ■ • , 5) of the order of unity, we require that x <C 1, y *C 1, and y 3> x in the asymptotic past. 
Again we shall set — in order to avoid that the last term in Eq. (28) always dominates the cosmological dynamics. 
At the dS point we have the following relations 

/j - ^2 24(2p 2 x 2 s + 5px ds + 9) 9 . 48 - 

e2(d 2 px dS ) = ^ — — + FJ a -TT^' ( 30 ) 

pxds + 9 64 3125 

(j - ,3 81(px ds + 2)(p£ ds +3) 243 , 243 - 

{d 3 px ds ) = — + —a - ——/3 , (31) 

pxds + 9 256 3125 



where 



e 4 (d 4 pi ds ) 4 , /3 = (d 5 px dS ) 5 ■ (32) 



In the above two theories we wish to study 1) whether different Galileon-like actions have some common feature, 
and 2) how they differ. 



IV. CONDITIONS FOR THE AVOIDANCE OF GHOSTS AND LAPLACIAN INSTABILITIES 

In order to discuss the stability of theories described by the Lagrangian (9) in the cosmological context, it is 
necessary to study linear perturbation theory on the flat FLRW background. Let us consider the perturbed metric 

ds 2 = -[1 + 2#(i, sc)] dt 2 + dix{t, x) dtdx 1 + a 2 (t)[l + 2$(<, x)] dx 2 , (33) 

where fy, <&, and \ are scalar metric perturbations. We have chosen the gauge 6(f) = without a non-diagonal 
scalar perturbation in the spatial part of the metric, i.e. dijj = (see Refs. [42] for the details of gauge-invariant 
cosmological perturbation theory). Taking into account two perfect fluids with the equations of state Wi = Pi/ pi 
(i = 1,2), there are three propagating scalar degrees of freedom. The velocity potentials Vi (i = 1, 2) are related with 

the energy-momentum tensor Tj 1 ^, as = —{p% + Pi)djVi (i = 1,2). 

Expanding the action (8) at second-order in perturbations, we find that the field ^ can be integrated out together 
with x- Introducing the vector Q = (yi,V2,$), we obtain the following second-order action for scalar perturbations 
(see Refs. [32, 43] for the details of such analysis): 



SS {2) = - I dtd 3 xa 3 



(34) 
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where A, C and D are 3x3 symmetric matrices and B is an antisymmetric matrix (we do not write explicit forms 
for them). 

Let us consider tensor perturbations with 6gij = a 2 hij, where hij is traceless (h l i = 0) and divergence-free j = 
0) . We also expand the action (8) at second-order in terms of the two polarization modes, hij = ft.® e® + h® efj , where 
e® and e® are the polarization tensors. For the polarization mode /i® , the perturbed action is given by 

,,2 



ssft = i / 



(35) 

where we will show the explicit forms of Qt and c T later. The conditions for the avoidance of ghosts and Laplacian 
instabilities of tensor perturbations correspond to Qt > and c T > 0, respectively. Note that the same expression 
also holds for h®. 

In Sees. IV A and IV B we study the general theories described by the Lagrangian (9) without imposing that /i 
and £ are zero. In Sec. IV C we shall apply our formula to specific theories with /i = and £ = 0. 

A. No-ghost conditions 

In order to avoid that a ghost mode appears for scalar perturbations, the matrix A needs to be positive definite. 
This leads to the following three no-ghost conditions 

pi > , (36) 

Wi 
1 + W2 

n . 7i 4(pi + p2 + 72)71 ~ 9t| ^ ( 

Qs = y -2 >0 ' ( 38 ) 



P2 > , (37) 



where 

71 = -6F - 9e 4 ./4 + l&Hfl - 24H£ , (39) 

72 = -9Fi7 2 - 9HF - j l /2 - 3e 2 / 2 2 /2 + 15H/| - 315e 4 ^ 2 /4/2 + 189i? 3 /| ~ 60H 3 ^ , (40) 

73 = - 2F + 2/| - 30e 4 i?/ 4 4 + 42i? 2 /| - 24i/ 2 ^ . (41) 

For radiation (w\ = 1/3) and non-relativistic matter (u>2 — +0) the conditions (36) and (37) are automatically 
satisfied. Hence we only need to consider the condition (38) to avoid the appearance of ghosts for scalar perturbations. 
The no-ghost condition for tensor perturbations is given by 

Qt = \f+ \e A ff - ^Hf! + 2Hi > . (42) 

B. Conditions for the avoidance of Laplacian instabilities 

The propagation speeds cs of three scalar degrees of freedom is known by solving the equation 

det(c|A - C) = . (43) 

The standard velocities for two perfect fluids correspond to c 2 s = W\ and Cj = W2, which are positive for both radiation 
and non-relativistic matter. The stability condition coming from the third solution is given by 



c 2 = 27 2 [73 - 2(1 + Wl )pi - 2(1 + w 2 )p2 - 73#] ~ 4717173 + 7 2 7 4 > Q 
C ' S 7i[4(pi +P2+ 72)71 - 973] 

where 



74 = - 18F + 9e 4 / 4 4 + 54/ 5 4 /5 - 72| . (45) 



7 



Under the no- ghost condition (38), this translates to 

2 7l 2 [73 - 2(1 + w x )pi - 2(1 + w 2 )p 2 - 73#] - 47i7i73 + 7a74 > . (46) 
The Laplacian instability of tensor perturbations is absent provided that the propagation speed squared is positive: 

4= 2Jr -^- 6 ^ + 8| .>0. (47) 
C. de Sitter stability 

On the dS background there are no matter fields, so that only one scalar mode propagates. The second-order action 
for scalar perturbations at the dS fixed point is 



Ss s ] = \ I 



* 2 - -§ (v$) 5 



a* 



(48) 



where Qs and c| correspond to those given in Eqs. (38) and (44) with the limits H — > H^s and px,2 — > 0. Let us 
discuss the conditions for the avoidance of ghosts and instabilities on the dS solutions we have discussed in Sec. III. 
We shall consider two theories described by the functions (17) and (26), with fx = and £ = 0. 

1. Constant F 

For the theory we discussed in Sec. Ill A, the conditions (38), (44), (42), and (47) reduce to 

Qs _ 4-9(a-2/3) 2 



M 2 " 3(a- 2/3)2 



>0, (49) 



4 = "^-97a-2fl»] > °' (50) 



(a - 2/3) (4 + 15a 2 - 48a/3 + 36/3 
|jT = ^(2 + 3a-6/3)>0, (51) 

In Sec. V we will show the allowed parameter space in the (a, /3) plane after deriving other conditions. 
From the action (48) we obtain the equation for $ in Fourier space: 

sbs (*>•*) +< *?*- - < 53 > 

where k is a comoving wavenumber. The solution for the homogeneous perturbation (k = 0) is 

^cx + czf-^at, (54) 

where c\ and C2 are integration constants. Since Qs is constant on the dS solution, the second term on the right hand 
side of Eq. (54) decays with time by noting that the scale factor evolves asaoc e Hdst . For the same reason, the tensor 
perturbation remains always stable in the limit k — > 0. This means that the dS fixed point is always classically stable 
under homogeneous perturbations. 

2, Non-constant F 

For the theory we discussed in Sec. IIIB, the conditions (42) and (52) reduce to 



i' r 



1 3 _ 3 - 
a - — — -8 



2 1024 6250 

4[15625<5 + 384/3(pi dS - 5); 
15625(3a + 512) - 7680/3 



>0, (55) 
>0. (56) 



s 



The expressions for Qs and c| are more involved, but Qs is proportional to </> p as in the case of Qt- Integrating the 
relation <j>/{H<j>) = x ds = constant, it follows that <j> oc a £dS . Since a 3 Q s oc e (3+pi d s)ffdst ) we g n( j f rom Eq (54) t h a t 
the homogenous perturbation $ evolves as 



$ = Si +5 2 e^ (3+p5ds)ffdS * , 
where ci and £2 are constants. Hence the dS point is classically stable for 

3 + pids > . 



(57) 



(58) 



The stability condition (58) is satisfied for |p£ds| <C 1. In this regime the conditions (38) and (44) for the scalar 
perturbation reduce to 



Qs 



M p i 

pxds 
27 



243 



+ O(x%)>0 



> 0. 



(59) 
(60) 



For positive <j> the no-ghost condition (59) is satisfied. If px^s < 0, the Laplacian instability of the scalar mode can 
be avoided. 



V. COSMOLOGY BASED ON THE COVARIANT GALILEON THEORY 

First we study cosmological dynamics for the covariant Galileon theory described by the functions (17) with = 0. 
This was partially discussed in the letter [34], but in this paper we shall thoroughly study the cosmology in such a 
theory with detailed numerical simulations. 

In the presence of radiation (pi = p r ,wi = 1/3) and non-relativistic matter (^2 = Pm,W2 — +0), we obtain the 
background equations from Eqs. (11) and (12): 

3M^H 2 = PDE + p r + p m , (61) 

Ml^H 2 + 2M 2 l H = -P DE - p r /3 , (62) 

where 

Pde = -e 2 dl^ 2 /2 + 3d 3 3 H<j) 3 /M 3 - A5e i diH 2 ^/(2M 6 ) + 21d 5 5 H 3 <p 5 /M 9 , (63) 

Pde = -e 2 d 2 2 4> 2 /2 - d 3 <j) 2 4>/M 3 + 3e 4 d^ 3 [8F0 + (3H 2 + 2H)<j>]/ (2M 6 ) - M 5 5 H^[5H4> + 2{H 2 + H)<j>]/M 9 . (64) 

The continuity equations for radiation and non-relativistic matter are given, respectively, by 

p r + 4Hp r = , p rn + 3H p m = . (65) 

From Eqs. (61), (62), and (65) the dark component also obeys the continuity equation 

PDB + 3H{pDE + PbE)=0. (66) 

We define the dark energy equation of state wde and the effective equation of state w c g, as 

wde= , W eff = -1 - rTTo , (67) 

PDE 3H Z 



where the latter is known by the background expansion history of the Universe. Using Eq. (66) together with the 
relation pde = 3Afp[iJ 2 S7DE, it follows that 

WDE = w cS - , (68) 

o"DE 

where a prime represents a derivative with respect to N = In a. 
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Each term in Eq. (20) has the difference of the order of x/y. The highest-order term in £1de comes from the term 
£5, i.e., of the order of y 3 /x. For the dynamical analysis given below, it is convenient to introduce the following 
quantities 



r i = — = — [-FT . r 2 = — 5- = — 3 • ( 69 ) 



At the dS fixed point one has n = 1 and r 2 = 1. In terms of n and r2 the density parameter (16) can be written as 

^de = ~(3a - 4/3 + 2)r 3 r 2 + (9a - 9/3 + 2)r 2 r 2 - ^-a ri r 2 + 70r 2 , (70) 

where a and /3 are defined in Eq. (25). Here we have employed the relations (23) and (24) to eliminate the terms e 2 d 2 
and d 3 . 

It is convenient to use the variables a and /3 for several reasons. First, the coefficients of physical quantities [such 
as 51de in Eq. (70)], autonomous equations, quantities related with no-ghost and stability conditions can be expressed 
in terms of a and ft. Second, the equations of motion, together with linear perturbation theory, are not subject to 
change under the following change of parameters x^s —> 7^dS an d di — > di/l (with i = 2,3,4,5), where 7 is a real 
number. In this case, depending on the parameter 7, there are infinite choices for the Lagrangian coefficients di that 
lead to the identical physics for the same a and ft. Therefore, constraining the parameter space in terms of a and 
ft allows us to remove the arbitrariness of the 7 rescaling. This also shows that one can set Xds = 1 without losing 
generalities. 

If n <C 1 at early times, the highest-order term £5 gives the dominant contribution to the dark energy density 
firjE- In this case it is expected that the cosmological Vainshtein mechanism can be at work to recover the General 
Relativistic behavior. If n 1 initially, the dominant contribution to firjE comes from the term C 2 . In this case the 
field energy density decreases rapidly as in the standard massless scalar field and hence the solutions do not approach 
the dS fixed point at late times. 

The conditions (38), (42), (44), and (47) for the avoidance of ghosts and instabilities of scalar and tensor perturba- 
tions reduce to 

Qs 6(1 + fJLl)((il + ji-i + ~ 2/x 3 - M3) > q / 71 n 



M p 2 ! (1 + M3) 2 

Qt 1 , 3 3 

(1 + Mi) 2 [2 M3 - (1 + fx 3 )(5 + 3w cS ) + 4fi r + m m ] - 4/4(1 + /n)(l + M2 ) + 2(1 + M3 ) 2 (l + M4 ) 



m -- + -a ri r 2 --ftr 2 >Q, (72) 



6(1 + ni)((ii + /i 2 + jUi/i 2 - 2^3 - 



> , (73) 



2 _ 2ri(2 - arir 2 ) - 3ft(r 2 r[ + nr 2 ) 



2r 1 (2 + 3ar 1 r 2 -6/3r 2 ) >U ' (74j 
where 

Hx = 3anr 2 /2 - 3ftr 2 , (75) 

fi 2 = (3a - 4/3 + 2)r 3 r 2 /2 - 2{9a - 9ft + 2)r\r 2 + A5anr 2 /2 - 28ftr 2 , (76) 

/i 3 = -(9a - 9ft + 2)r 2 r 2 /2 + 15ar!r 2 /2 - 21/3r 2 /2 , (77) 

IH = -anr 2 /2 - 3ftr 2 (r[/r 1 + r 2 /r 2 )/4 . (78) 

From Eqs. (11), (12), and (13) we obtain the following differential equations for the variables n, r 2 , and f2 r : 

rj = ^■(r 1 -l)r 1 [r 1 (r 1 (-3a + 4/3-2) + 6a-5/3)-5/3] 

x [2 (O r + 9) + 3r 2 (r?(-3a + 4/3 - 2) + 2r 2 (9a - 9/3 + 2) - 15na + 14/3)] , (79) 

r 2 = -^[r 2 (6r 2 (r 2 (45a 2 - 4(9a + 2)/? + 36/3 2 ) - (Q r - 7)(9a - 9/3 + 2)) + r?(-2(Q r + 33)(3a -4/3 + 2) 

-3r 2 (-2(201a + 89)/? + 15a(9a + 2) + 356/3 2 )) - 3ria(-28fi r + 123r 2 /3 + 36) + 10/3(-llO r + 2lr 2 ft - 3) 
+3r^r 2 (9a 2 - 30a(4ft + 1) + 2(2 - 9ft) 2 ) + 3r\r 2 (3a - Aft + 2) 2 + 3rjY 2 (9a - 9/3 + 2)(3a -4/3+ 2))], (80) 

Q' r = -|a.[r 2 (4(ft r - l)(9a - 9ft + 2) + 6r 2 (-15a 2 + 36a/3 + 4(2 - 9ft)ft)) - 2r?((fi r - l)(3a -4/3 + 2) 

+9r 2 (18(a + 1)0 + a(9a + 2) - 36ft 2 )) + 12na(-3fi r + 22r 2 /3 + 3) - 10/3(-4f2 r + 21r 2 /3 + 4) 

+rj t r 2 (549a 2 + a(330 - 840/3) + 2(2 - 9/3) 2 ) + 3rfr 2 (3a - 4/3 + 2) 2 - 12rjY 2 (9a - 9/3 + 2) (3a - 4/3 + 2)], (81) 
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where 

A = 2r*r 2 [72a 2 + 30a(l - 5/3) + (2 - 9/3) 2 ] + 4r 2 [9r 2 (5a 2 + 9a0 + (2 - 9/?)/3) + 2(9a - 9/3 + 2)] 

+4r?[-3r 2 (-2(15a + 1)0 + 3a(9a + 2) + 4/3 2 ) - 3a + 4^ - 2] - 24r ia (16r 2 /3 + 3) + lO0(21r 2 + 8) . (82) 

The Hubble parameter obeys the following equation 

H' _ 5r[ r' 2 

where r'i/ri and r 2 /r 2 are known from Eqs. (79) and (80). 

A. Tracker solutions (n = 1) 

From Eq. (79) we find that there is an equilibrium point characterized by 

n = l , (84) 

at which the density parameter (70) reduces to 

Ode - r 2 . (85) 

From Eq. (69) we find that xH 2 = constant along the solution (84). Hence the field velocity evolves as 

<j) oc H' 1 , (86) 

which has the dependence <f> oc t during the radiation and matter eras. Since the field is effectively frozen at early 
times, this shows the implementation of the cosmological Vainshtein mechanism. 
Along the solution (84). the other two equations can be written as follows 

2r 2 (3 - 3r 2 + n r ) 

^ = 1 + T2 ' (87) 

= »r ("r ~ 1 ~ 7r 2 ) 

1 + r 2 

which do not depend on a and 0. We then have the following three fixed points 

(A) (ri,r2,fi r ) = (1,0,1), (B) (r u r 2 ,n r ) = (1,0,0), (C) (n, r 2 , O r ) = (1, 1, 0) . (89) 

The points (A) and (B) can be realized during the radiation and matter eras, respectively, whereas the point (C) 
corresponds to the dS solution. 

The stabilities of these fixed points can be analyzed by considering linear perturbations Sri, Sr 2l and 5fl r about 
them. For example, the perturbation Sri satisfies 

, 9 + q. + 3r 2 
frl = ~ 2(1 +r 2 ) (90) 

This shows that, in the regime < r 2 < 1 and fi r > 0, the solution is stable in the direction of r\. Defining the vector 
Sr = t (Sri, 5r 2 , Srs), one can write the perturbation equations in the form 

Sr' = MSr, (91) 

where M. is the 3x3 matrix. The eigenvalues of the matrix Ai for the points (A), (B), (C) are given by 

(A) (8, 1, -5) , (B) (6,-1, -9/2) , (C) (-3, -3, -4) . (92) 

This shows that (A) and (B) are saddle, while (C) is stable. Hence the solutions finally approach the stable dS point 
(C). This dS stability is consistent with the analysis in Sec. IV C 1 based on homogeneous perturbations. The solution 
(84) can be regarded as a tracker that attracts solutions with different initial conditions to a common trajectory. 
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Along the tracker we have p D E = 3M 6 /H 2 , P DE = -3M 6 {2 + w eS )/H 2 , and 

n r + 6 n r - 6r 2 

wde = -2 - w cS = -— — — , w eS = — — . (93) 

3(r 2 + 1) 3(r 2 + 1) 

During the cosmological sequence of radiation, matter, and dS eras the dark energy equation of state evolves as 
u>de = —7/3 — > —2 — > —1, whereas the evolution of the effective equation of state is w c s = 1/3— > — > — 1. This 
peculiar evolution of u>de can be useful to constrain the covariant Galileon theory from observations. 

Equations (87) and (88) are simple enough to be solved analytically. In fact, combining Eqs. (87) and (88), it 
follows that 



2 



2^ , (94) 



7*2 fi 

which has the solution 

r 2 = Cl a 8 n 2 r , (95) 

where c\ is a constant of integration. Substituting this solution into Eq. (88), we find two branches that differ from 
each other in the early cosmological limit. The viable branch of solutions is given by 



c 2 a - 1 + VI - 2c 2 a + c 2 2 a 2 + 4 Cl a s 
^ = 2^ ' (%) 

where c 2 is another constant. Since fl r ~ 1 + c 2 a at early times (a <C 1), we require that c 2 < (provided Ode > 0). 

The coefficients Ci and c 2 can be found by using the present density parameters of radiation and non-relativistic 
matter, i.e. O r (a = 1) = Of 5 and O m (a = 1) = £1$ . Using the relation (85) as well, we find 

1 - o (0) - o (0) o (0) 

c l = m > c 2 = TTvT ■ I 97 ) 

The density parameter of dark energy evolves as 

Ode = c ia 8 n 2 r . (98) 

Hence the density parameters Ode, O n and O m as well as ijjde and u; e g are analytically known in terms of the 
function of a (or the redshift z = l/o— 1). 

At the dS point (C) the conditions for the avoidance of ghosts and instabilities have been already estimated in 
Eqs. (49)-(52). Let us consider the points (A) and (B), which are characterized by r\ = 1 and r 2 <C 1. In this case 
Eqs. (71)-(74) are simplified to give 

Q s /Mp 1 -3(2-3a + 6/3)r 2 >0, (99) 
Qr/Mpi = 1/2 + 3(o - 2/3)r 2 /4 > , (100) 
2 8 + lQ a -9/3 + O r (2 + 3a-3/3) 

3(2 -3a + 6/3) >U ' (101) 

4 ~ 1 - (4a + 3^ + 3/3 r 2 r )r 2 /2 > 0. (102) 

Since r 2 <C 1 the conditions (100) and (102) are automatically satisfied. From Eq. (99) the sign change of r 2 means 
the appearance of the scalar ghost. If we choose the initial conditions with r 2 > 0, then Eq. (99) requires that 

2 - 3a + 6/3 > . (103) 

Let us consider the intermediate regime between r 2 -c 1 and r 2 = 1. As long as the conditions Qs > 0, Qt > 0, 
and c| > are satisfied both in the regimes r 2 <C 1 and r 2 = 1, the violation of these conditions does not occur in 
the intermediate epoch. However, the tensor propagation speed squared can be negative even if the conditions (52) 
and (102) are satisfied. Along the tracker Eq. (74) gives 

2 _ 2 + (2 - a - 9/3 - 3/3Q r )r 2 + (9/3 - a)r 2 

(l + r 2 )[2 + 3(a-2/3)r 2 ] ' (HMJ 
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FIG. 1: The viable parameter space in the {a, 13) plane determined by the conditions (49)-(52), (101), (103), and (106) along 
the tracker solution r\ = 1 [34]. 

The transition to the dS solution (r 2 = 1) occurs only recently, so that the term £l r can be neglected in Eq. (104). 
Then c 2 ^ has an extremum at 

_ 4a - 15/3 ± 3a//3(30/3 - 8a + 12a 2 - lba(3 - 18/3 2 ) 
V2 ~ 15/3 - 4a + 27a{3 - 54^ 2 ' (1 ' 

If a = 1.9 and /3 = 0.8, for example, the physical solution corresponds to the plus sign in Eq. (105), i.e. r2 = 0.636, 
at which c 2 ^ has a minimum. As /3 approaches 1, the minimum values of c 2 ^ get smaller. For /3 around 1, c 2 ^ can be 
negative for the plus sign of Eq. (105). This leads to the following condition for f3 > 0: 

2/3 < a < 12^/3- 9/3- 2. (106) 

If (3 < 0, then cf, remains to be positive. Hence we do not have any additional constraint in the regime f3 < 0. 

In Fig. 1 we plot the parameter space constrained by the conditions (49)-(52), (101), (103), and (106). For the 
solutions that start from initial conditions with n l,r2 *C 1 and then approach the dS attractor with r\ = 1 and 
r<i = 1, the parameters a and /3 need to be inside the purple region in Fig. 1. There is another case in which both r\ 
and r2 are initially much smaller than 1. We shall address this case in the next subsection. 



B. Solutions driven by the term £5 (n <l,rj< 1) 

From Eq. (79), it is clear that another equilibrium point exists, namely, r\ = 0. Let us now discuss this equilibrium 
point in more detail. In this case Eqs. (80) and (81) reduce to 

_ r 2 (21r 2 /j-liq.-3) 

21r 2 /3 + 8 ' (1U7j 

2q.(21r 2 /3-4a. + 4) rinR x 

a " 2i7^T8 ' (108) 

which depend on /3. The dominant contribution to the field energy density comes from the term £5, i.e. Ode = 7/3?"2- 
We have the following fixed points 

(A') (ri,ra,n r ) = (0,0 > l), (B') (ri,r 2 ,O r ) = (0,0,0) , (C) (r l5 r 2 ,Q r ) = (0, 1/(7/3), 0) , (109) 
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which represent radiation, matter, and dark energy dominated points, respectively. Perturbing Eq. (79) on the n = 
solution leads to 

21r 2 /3 + Q r + 9 

frl= 21r 2 /3 + 8 (110) 

which implies that none of the fixed points (A')-(C') can be stable. In particular the eigenvalues of the matrix A4, 
where Sr' — AiSr and Sr — t {Sr-y,8r2,Sr^), are given by 

(A') (5/4, 7/4, 1) , (B') (9/8, 3/8,-1) , (C) (12/11, -3/11, -14/11) . (Ill) 

This shows that the point (A') is unstable, whereas the other two are saddle. Recalling that the dS fixed point (C) 
discussed in the previous subsection is stable against homogenous perturbations, the solutions finally approach (C) 
instead of (C). Unless n is initially very small such that the solutions reach n = 1 only at late times, the system 
approaches the stable n = 1 direction much before the dS epoch. 

In the regime r\ <C 1 and r 2 -C 1 it is possible to derive analytic solutions for T\ and r 2 as well as for wde and w e g. 
In fact, Eqs. (79), (80), and (81) can be simplified as 

r , 1 ~i(fi r + 9)n > (112) 
r 2 ~i(im r + 3)r 2 , (113) 

n' T ~ -n r (i - n r ) , (114) 

where we have assumed that \(3\ is not very much smaller than unity. During the radiation domination (Q r = 1), 
integration of Eqs. (112) and (113) gives 

ricxa 5 / 4 , r 2 «a 7/4 , (115) 

whereas during the matter era one has 

ri oc a 9 / 8 , r 2 oc a 3 / 8 . (116) 

Eventually the solutions approach the tracker r\ = 1. 
In the regime ri <C 1, r 2 -C 1 one has 

WT)E ~ -(l + fi r )/8, w cff ~fi r /3. (117) 

This gives »de — — 1/4 and w c ff ~ 1/3 during the radiation era, whereas wde — — 1/8 and iy e ff — during the matter 
era. 

The condition (71) reduces to 

Qs/M^ ~ 60/3r 2 > . (118) 

The sign change of r 2 implies the appearance of ghosts. For the initial conditions with r 2 > we require that 

/8>0. (119) 

If the solutions start from the regime n <C 1, t% <C 1 and subsequently enter the regime n = 1, the allowed parameter 
space in Fig. 1 is restricted be /3 > 0. Since Qt/M^ ~ 1/2, the no-ghost condition for the tensor mode is automatically 
satisfied. 

The propagation speeds of scalar and tensor perturbations are given, respectively, by 

c|~ (l + fi r )/40, (120) 
c| ^ 1 + 3^2(5-30^/8, (121) 

which are both positive for < Q r < 1. The scalar mode remains sub-luminal during the radiation era (c 2 s = 1/20) 
and the matter era (c| = 1/40). Under the no-ghost condition (118) the tensor mode becomes super-luminal (although 
4 is very close to 1 ) . 
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FIG. 2: Evolution of Que, fi m , fir, and w e s versus the redshift z = 1/a — 1 for a = 0.3, /3 = 0.14, £2 = 1, £4 = 1, and x'ds = 1. 
We choose the initial conditions ri = 1.500 x 10~ 10 , r 2 = 2.667 x 10~ 12 , and fi r = 0.999992 at z = 3.63 x 10 8 . 

C. Numerical simulations for the cosmological dynamics 

Numerically we integrate Eqs. (79)-(81) to confirm the analytic estimation in the previous subsections. 

Let us consider the case in which the variables r\ and T2 are much smaller than 1 at the initial stage of cosmological 
evolution. Our numerical simulations show that r\ and ri evolve as Eq. (115) during the radiation era, whereas their 
evolution during the matter era is given by Eq. (116). Depending on the initial conditions of T\ and ti, the epoch at 
which the solutions approach the tracker (n = 1) is different. As we increase the initial ratio ri/r%, this epoch tends 
to occur earlier. After the solutions reach the tracker, the evolution of ri, O r , and Ode is given by Eqs. (95), (96), 
and (98), respectively. 

In Fig. 2 we plot one example for the evolution of density parameters Ode, O m , and O r as well as the effective 
equation of state w e f{. In this case the transition to the regime n — 1 occurs only recently, e.g., r\ = 0.99 around 
z = 0.07 with r2 ~ 0.6. After passing the present epoch, the solutions are attracted by the dS solution characterized 
by (ti,T2) = (1, 1)- Figure 2 shows that the sequence of radiation (O r = 1, w e s = 1/3), matter (O m = 1, w e s = 0), 
and dS (Ode = 1, w c s = —1) epochs is in fact realized. Unlike dark energy models based on f(R) theories, the 
Galileon model is not plagued by the presence of a rapidly oscillating mode associated with a heavy field mass in the 
early Universe. 

Figure 3 illustrates the variation of u>de for several different initial conditions and model parameters. The cases 
(a)-(d) correspond to a — 0.3, j3 — 0.14, £2 = 1, £4 = 1, and x^s = 1 with different initial conditions satisfying n < 1 
and f2 <C 1, whereas the case (e) shows the tracker solution starting from the initial condition ri = 1 and r-i <C 1 with 
the model parameters a = —1.5, (3 = —0.9, £2 = 1, £4 = — 1 5 and 2-ds = 1- Clearly the solutions with different initial 
conditions converge to the tracker, depending on the epoch at which the variable n grows to the order of 1. In the 
cases (a)-(d) the dark energy equation of state evolves as Eq. (117) in the regime r\ <C 1 and <§C 1 (»de — —1/4 
and wde — —1/8 during the radiation and matter eras, respectively), which is followed by the evolution given in 
Eq. (93) after the solutions reach the tracker at r\ = 1. As long as the tracking behavior occurs by today, the dark 
energy equation of state crosses the cosmological constant boundary (wde = — 1)- 

Numerically we find that for the initial conditions with r\ < 2 the solutions are typically attracted by the tracker. 
On the other hand, if n > 2, the system tends to approach the matter-dominated epoch with the growth of r±. In 
the latter case the dominant contribution to Ode comes from the term £2, so that Ode decreases as in quintessence 
without a potential. 

In Fig. 4 we plot the evolution of c 2 s for the same model parameters and initial conditions as those presented in 
Fig. 3. In the regime r\ <C 1 and r2 < 1, our numerical simulations in the cases (a)-(d) agree with the analytic 
estimation of the scalar propagation speed given in Eq. (120), i.e. c| ~ 1/20 and c| ~ 1/40 during the radiation and 
matter eras respectively. As the solutions reach the regime r\ ~ 1 with r% <C 1, c| approaches the value estimated by 
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FIG. 3: Variation of wue versus z for a = 0.3, /3 = 0.14, t-z = 1, 64 = 1, and Xds = 1 [cases (a)-(d)]. We choose four different 
initial conditions: (a) n = 5.000 x 10 -11 , r 2 = 8.000 x 10~ 12 , and 9, r = 0.999995 at z = 5.89 x 10 8 , (b) n = 1.500 x 10~ 10 , 
r 2 = 2.667 x 10" 12 , and Q. r = 0.999992 at z = 3.63 x 10 s , (c) n = 5.000 x 10" 9 , r 2 = 8.000 x 10~ 14 , and ft r = 0.99995 at 
2 = 6.72 x 10 7 , (d) n = 5.000 x 10" 6 , r 2 = 8.000 x 10~ 17 , and Q. r = 0.9986 at z = 2.04 x 10 6 . The case (e) corresponds 
to a = —1.5, /3 = —0.9, £2 = 1, £4 = —1, and Xds = 1 with initial conditions ri = 1, T2 = 10 -60 , and Sl r = 0.99999 at 
z = 3.12 x 10 s . 
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FIG. 4: Evolution of c| versus z for the same model parameters and initial conditions as given in Fig. 3. 



Eq. (101). When a = 0.3 and (3 = 0.14 the analytic estimation gives c| ~ 1.67 during the matter dominance, which 
agrees with the value at the plateau in the case (d) of Fig. 4. Finally the solutions reach the dS fixed point, at which 
c| shifts to the value given in Eq. (50), e.g.. c| = 1.01 x 10 -2 for a = 0.3 and /? = 0.14. 

For positive /3 one can show that under the conditions (49), (50), and (99) the scalar propagation speed estimated 
by Eq. (101) becomes super-luminal. However, the scalar mode can remain sub-luminal provided the solutions reach 
the regime r\ = 1 in the recent past. The cases (a) and (b) in Fig. 4 correspond to such examples in which the peak 
value of c| is smaller than 1. 
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FIG. 5: Evolution of versus z for two cases: (a) a = 0.3, /3 = 0.14, e 2 = 1, £4 = 1, and Xds = 1 with initial conditions 
n = 1.500 x 10~ 10 , r 2 = 2.667 x 10~ 12 , and Q, r = 0.999992 at z = 3.63 x 10 s , (b) a = -1.5, /3 = -0.9, e 2 = 1, e 4 = -1, and 
x ds = 1 with initial conditions r\ = 1, r 2 = 10~ 60 , and fi r = 0.99999 at z = 3.12 x 10 8 , and (c) a = 1.9, /3 = 0.8, e 2 = 1, 
e 4 = 1, and x dS = 1 with initial conditions n = 10~ 5 , r 2 = 10" 35 , and f2 r = 0.99999 at z = 3.12 x 10 8 . 

If /3 < there is a parameter space in which the scalar propagation speed (101) is sub-luminal, while satisfying the 
conditions (49), (50), and (99). In this case the initial conditions of t\ need to be close to 1. If t\ is smaller than 
the order of unity, the scalar ghost appears for negative /3. On the other hand, if T\ > 2, the solutions do not finally 
approach the dS fixed point. The case (e) in Fig. 4 corresponds to an example of the sub-luminal evolution of c 2 s for 
negative j3 with the initial condition r\ — 1. Since the solution stays on the tracker, the scalar propagation speed is 
given by Eq. (101) during the radiation and matter eras and by Eq. (50) at the dS point. 

For the initial conditions with ri < 1 and T2 <C 1 the tensor propagation speed starts to evolve from the value 
estimated by Eq. (121), which is slightly super- luminal under the no-ghost condition f3 > for the scalar mode. After 
the solutions reach the regime r\ ~ 1 and <C 1, is still close to 1 because it is described by Eq. (102). The tensor 
propagation speed finally approaches the value (52) at the dS point. During the transition from the regime r-i <C 1 to 
the regime r<i ~ 1, there is an epoch at which can have either the maximum or the minimum. In the case (a) of 
Fig. 5, the analytic formulas in Eqs. (104) and (105) show that has a minimum value 0.799 at r-i = 0.662 [plus sign 
of Eq. (105)], whereas in the case (b) c\ possesses a maximum value 1.690 at T2 — 0.412 [minus sign of Eq. (105)]. 
This estimation agrees well with the numerical results shown in Fig. 5. In the case (c) of Fig. 5 the condition (106) 
is violated, so that has a negative minimum. In the region where a and (3 are positive, the condition (106) needs 
to be satisfied to avoid the temporal Laplacian instability of the tensor mode. 

If the solutions start from the regime r± ~ 1 and T2 <C 1, then the tensor propagation speed (102) can be sub-luminal 
under the condition 4a + 3/3 + 3/3Q r > for the branch Ti > 0. In this case, however, exceeds 1 at the dS point, 
as long as the conditions (49)-(51) and (99)-(101) are satisfied. Since of* > 1 in the regime r% <C 1 and <C 1 as 
well, it is not possible to avoid the appearance of the super-luminal mode for tensor perturbations. However, the 
super-luminal propagation does not necessarily imply the inconsistency of Galileon theory because of the possibility 
for the absence of the closed causal curve [30] . 

VI. COSMOLOGY BASED ON THE MODELS WITH NON-CONSTANT FUNCTIONS F(<j>) 

We shall proceed to the cosmology for the theories with non-constant F in which the functions F and /j (i = 
1, 2, • ■ ■ ,5) are given in Eq. (26) with = 0. We take into account radiation (p\ = p r , wi = 1/3) and non-relativistic 
matter (p2 = Pm, ^2 = 0), which satisfy the continuity equations (65). Taking the time-derivative of Eq. (11) and 
combining it with Eq. (12), we obtain the equations of motion for <fi and H. Then the dimensionless variables x, y 
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defined in Eq. (29) and the radiation density parameter fJ r = p r /(SFH 2 ) obey the following equations 

x = £[29160000000000 x 3 (-8Q r (12 + d 2 2 e 2 px) + (4 + (-4 + d%e 2 )px)(24 + px(24 + d%e 2 px))) 
-1440000000000 d^p 2 x 3 (-216 - 216 fl r + px{24 + pi(-192 + d 2 2 e 2 (45 + pZ))))y 2 

-625000000 p 3 x 2 (2048 c^p 2 x 2 (-9 + pi) + 6561 d\e A (144 + 80 fl r + px(20 - px(-76 + dje 2 (19 + px)))))y A 
+583200000 p 4 £(3125 dld\e A p 2 x 2 ! (-ll + pi) + 192 df(600 + 280 O r + px(88 + px(272 
-d 2 2 e 2 (73 + 5px)))))i/ 6 - 84375 p 7 i 2 (11390625 d^(-6 +px) + 524288 (-45 + 2px))y 8 
+15746400000 d^e 4P 8 i(-59 + 9px)y 10 - 6019743744 4V(-5 + px)y 12 ]/(60A) , (122) 
y' = xy[116640000000000x 2 (24 - 240 r + px(-48 - 24pi + t^e 2 (12 + 5px)) - 4320000000000ci^ 2 £ 2 (72 - 24Vl r 
+px(-16 + (-32 + U\t 2 )px))y 2 - 625000000p 3 i(2048^p 3 i 3 + 2187d|e 4 (96fi r + pi(108 -px(-120 
+^e 2 (18 +pi)))))y 4 + 64800000p 4 (3125d 3 i c£[e 4 p 2 z 2 (-9 + lOpx) + 1728^(120 + 120fi r 
+pi(128 + p£(144 - d\t 2 (21 + 2pi)))))f - 253125p 7 £(3796875d 4 i e 4 : pi + 524288df^(-3 + px))f 
+5248800000^ e 4 p 8 (-9 + 28pi)y 10 - 60197437444°p 10 y 12 ]/(60A) , (123) 
n' r = pft r [972000000000x 4 (-48pi + ^e 2 (8 - 8ft r +px(16 + (-4 + d%e 2 )px))) - 48000000000dfpx 3 (144 - 144fi r 
+px(-72 + px(-l68 + dle 2 (36 + px))))y 2 + 62500000p 2 i^(-1024c^p 2 x 2 (-2 +px) + 729d|e 4 (144 - 144fi r 
+px (-192 + p£(-204 + ^e 2 (43 + 3p£)))))y 4 + 2160000p 3 x(3125d? ; d 4 : e 4 p 2 i 2 (-18 + 17px ) + 1728d|(-160 
+160ft r +px(280 +p£(248 - t^e 2 (54 + 5px)))))f - 5625p 6 i 2 (11390625dfe 4 ! (l + px) 

+65536d 3 ^(-36 + 23px))y 8 + 174960000d 4 dfe 4 p 7 x(88 + 63px)y 10 - 5016453124°p 8 (2 +px)y 12 ]/A , (124) 

where 

A = p[5016453124°p 8 y 12 - 3732480000dfp 3 xy 6 (-160 + M 4 e 4 p 4 y 4 ) - 62500000x 4 (-1024(-27 + d\p 2 y 2 ) 2 

+243^62(512 + 3dje iP 4 y 4 )) + 151875p 2 f 2 y 4 (65536d|p 2 j/ 2 (-45 + d\p 2 y 2 ) + 84375^e 4 (-512 + 5d|e 4 p 4 y 4 )) 
+4320000pi 3 y 2 (108p 2 y 2 (9375ci|e 4 + \M\d\t 2 p 2 y 2 ) - 3125d|(-512 + 9^e 4 p 4 y 4 ))] . (125) 

The dS fixed point with x = Xds = y and Cl r = exists under the conditions (30) and (31). Since the theory has 
a nonminimal coupling F(<p)R, it is possible to place constraints on the values of x around today from the variation 
of the effective gravitational coupling, G e g = [87rF(</))] _1 . The Lunar Laser Ranging experiments give the bound 
|Goff/Gcff| < 1-3 x 10~ 12 yr _1 [44], or in terms of the present Hubble parameter H , \G e s/G e g\ < 0.02if [45]. In our 
theory \G e e/G e s\ = \px\ H, which gives the constraint |px| < 0.02 around today. Since the value of ids is not much 
different from x today, we employ the following criterion 

|px dS | < 0(0.01) . (126) 

Under this bound the condition (58) is always satisfied, which means that the dS solution is classically stable. From 
Eq. (60) the Laplacian instability of scalar perturbations at the dS point can be avoided for pxds < 0. The no-ghost 
condition (59) is satisfied provided that {<p/M p \) p > 0. 



A. Initial conditions with y 2 S> |x| 

If y 2 3> |x| in the early cosmological epoch, then the term £5 dominates over the terms £2,3,4, i.e. ^de — 
—pi + (7/3125)dfp 5 j/ 6 /x. In order to avoid the dominance of dark energy during the radiation and matter eras we 
require that |px| <C 1 and |d|p 5 {/ 6 | <C |x|. In this regime the quantities Qs and Qt defined in Eqs. (38) and (42) are 
approximately given by 

Qs „ 12 d 5 5 p 5 y 6 f * V Q T 1 / \ p f 3 4^f \ 

M 2 625 x \M P J ' M 2 2 \M P J \ 3125 x J ' { ' 

The tensor ghost is absent for {(j)/M p \) p > 0. Since the evolution of the field is given by <f> = (pi exp(f^ xdN), where 
(pi is the initial field value at N = Ni, the condition ((p/M p \) p > is satisfied for </>,; > 0. For the avoidance of the 
scalar ghost we require that 



d^px > . 



(128) 
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FIG. 6: Evolution of Ude, fim, ^r, and w e ff versus the redshift z for the model with p = 1, 62 = 1, £4 = — 1, ^4 = 1, <fe = 1, 
= 0, and ids = 0.007. The initial conditions are chosen to be x — 1.0 x 10~ 18 , y = 1.5 X 10~ 5 , and Sl r = 0.99992 at 
2 = 3.9x 10 7 . 



In the regime y 2 3> |5| and |c?|p 5 y 6 | <C |i| the scalar and tensor propagation speeds defined in Eqs. (44) and (47) can 
be estimated as 



4 



1 



3 „ _ ^l/y 6 



25000 



(4-3fi r )- 



(129) 



Since dp is close to 1, the tensor instability can be avoided. If d§ > 0, then there is no instability for the scalar 
perturbation (c| > 0). In the regime x 2 -C | ^|j3 4 y 6 1 we have c| ~ (1 + il r )/40 > 0. If x 2 > \d^p 4 y 6 \, it can happen 
that the scalar perturbation is subject to the Laplacian instability for negative dg. 

In the regime y 2 ^> \x\ and |c?|p 5 y 6 | -C |S| the autonomous equations (122)-(124) are simplified as 



y' 



8' 



15 + 70 r + 625(1 - n r 
625 

l + O r + — (l-OO 



3 



-fi r (l - Or) . 



(130) 

(131) 
(132) 



From Eq. (132) there are two fixed points characterized by O r = 1 and Q r = 0. As long as the condition x 2 <C |d|p 4 ?/ 6 | 
is satisfied, the evolution of the variables x and y during the radiation era (f2 r = 1) is given by 



11/4 



y (x a, 



3/4 



(133) 



whereas during the matter era (Q r = 0) one has 



x oc a 



15/8 



y oc a 



3/8 



(134) 



In both cases x grows faster than y. If the quantity x 2 / (d^p^y 6 ) becomes larger than the order of unity, the evolution 
of x and y is subject to change. 

For the solutions starting from the regime \y 2 \ 3> \x\ the condition (128) needs to be satisfied initially. Then there 
are two possible cases: (i) pi > and d$ > 0, and (ii) pi < and d$ < 0. The avoidance of the scalar Laplacian 
instability at the future dS fixed point requires that px^s < 0. However, if we demand the viable cosmology by today 
(the redshift z > 0), the condition px^s < is not necessarily mandatory. In general, if the variable x changes its 
sign during the cosmic expansion history, this signals the violation of the conditions for no ghosts and no Laplacian 
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FIG. 7: Variation of x, y, and Q r versus z for the same model parameters and the initial conditions as those given in Fig. 6. 
The solution finally approaches the de Sitter attractor with x = x^s = 0.007. 




log 1Q (1+Z) 

FIG. 8: Evolution of the dimensionless variables Qs = 2Q s / (M^~ p 4> p ) and Q T = 2Q T /{M*~ P 4> p ) versus z for the same model 
parameters and the initial conditions as those given in Fig. 6. The signs of Qs and Qt remain to be positive. 

instabilities. For example, this can be seen in the expression of Qs and c| in Eqs. (127) and (129) in the past 
asymptotic regime. In fact we have numerically confirmed the violation of at least one of those conditions. In the 
following we shall study the cosmological dynamics in which the sign of px at the early epoch is same as that of pi ds ■ 
In the case (i) the condition pxds < is violated, but it is possible to realize cosmological trajectories in which all 
the required conditions are satisfied by today. In the case (ii) the condition pids < is met, but we need to check 
whether there are no violations of the no-ghost and stability conditions in the cosmic expansion history. 

Let us first discuss the cosmological dynamics in the case (i) with px^s > 0. In Fig. 6 we plot the variation of Ode, 
Q m , Cl r , and w e s for the model with p = 1, 62 = 1, £4 = — 1, d± = 1, ds = 1, = 0, and x^s = 0.007 (in which case 
the condition (126) is satisfied). The constants cfe and ^3 are known from Eqs. (30) and (31). We choose the initial 
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FIG. 9: (Left) Evolution of c| and Sp versus z for the same model parameters and the initial conditions as those given in 
Fig. 6. (Right) The enlarged version for the evolution of c| in the regime —1.5 < log 10 (l + z) < 0. The field propagation speed 
becomes negative in the future (around z m —0.87). 



conditions x = 1.0 x 10~ 18 , y = 1.5 X 10~ 5 , and £l r = 0.99992 at the redshift z = 3.9 x 10 7 , in which case \y 2 \ > \x\ 
and |dgp 5 y 6 | <C |i| initially. The background evolution in Fig. 6 shows that the sequence of radiation, matter, and dS 
eras is realized in this case. 

Figure 7 illustrates the evolution of the variables x and y as well as f2 r . We find that i approaches the dS attractor 
with ids = 0.007 without changing its sign. In the regime i 2 <C |c?5P 4 y 6 | the evolution of i and y is well described 
by the analytic estimation (133) during the radiation era. However, around z < 10 5 , the last terms in Eqs. (130) and 
(131) starts to give rise to the contribution to the evolution of i and y. As we see in Fig. 7, i and y evolve differently 
from the analytic estimation (133) and (134) for z < 10 5 . 

In Figs. 8 and 9 we plot the variation of the quantities Qs = 2Qs / (M^ p (j) p ), Qt = %Qt / '(M^ p '(jf), c|, and c T for 
the same model parameters and initial conditions as those given in Fig. 6. We find that Qs grows rapidly, whereas 
Qt is always close to 1. Since both Qs and Qt are positive, the appearance of the scalar and tensor ghosts is avoided 
in this case. 

Figure 9 shows that c| starts to evolve from the value around 0.05, as estimated analytically in Eq. (129). For 
z < 10 5 the contribution of the second term in the expression of c| in Eq. (129) becomes important, which leads to 
the increase of c|. For the model parameters given in Fig. 6 the scalar propagation speed slightly exceeds 1 during 
the transition from the matter era to the dS epoch. In Fig. 9 we find that c| remains positive until recently (z > 0). 
However, since the sign of i is always positive, c| is negative at the dS point, i.e. c| ~ — pids/27 = —2.6 x 10 -4 . 
The crossing of c s at occurs in future around the redshift z sa —0.87. The tensor propagation speed squared is 
always close to 1 (slightly larger than 1), which means that the Laplacian instability of the tensor perturbation can 
be avoided. 

Let us next discuss the case with (ii), i.e. px < and d$ < initially. In Fig. 10 we plot one example for the 
evolution of c| with p = 1, = —1, and ids = —0.007. In this case the density parameters as well as the effective 
equation of state evolve similarly as those in Fig. 6. However, even if the variable x starts from negative values, x 
crosses for many times before reaching the dS solution with ids = —0.007. As we see in Fig. 10, this leads to the 
violation of the condition c 2 s > by today. In addition the quantity Qs also becomes negative during some periods. 
We have run our numerical code for many other cases in which the condition pxds < is satisfied and found that 
in the case (ii) it is difficult to find a viable cosmological trajectory along which all of the no-ghost and stability 
conditions are satisfied. 

In summary, we have shown that the cosmological solutions along which px > and d$ > initially and pids > at 
the de Sitter attractor can evade the problems of the ghost and instability problems for z > 0. In this case, although 
the scalar Laplacian instability is present at the de Sitter fixed point, the crossing of c 2 s at occurs at some time in 
future. We have also run the numerical code for the initial conditions with \y 2 \ > |i| and found similar properties of 
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FIG. 10: Evolution of c| versus z for the model with p = 1, es = 1, d = — 1, cLt = 1, ds = — 1, d{ = 0, and ids = —0.007. The 
initial conditions are chosen to be x = -4.5 x 10 -15 , y = -7.0 x 10~ 5 , and Q, r = 0.996 at z = 1.56 x 10 6 . In this case the scalar 
perturbation is subject to the Laplacian instability for many times by today. 

solutions to those discussed in this section. 

B. Initial conditions with y 2 \x\ 

Finally we shall study the case in which \y 2 \ <C \x\ in the early cosmological epoch. In this regime the term £2 is 
the dominant contribution to J7de relative to £3,4,5, i.e. Qde — — px — e 2 d 2 P 2 x 2 /24. The quantities Qs and Qt are 
approximately given by 

Qs _ p 2 x 2 (6-e 2 4) ( <t> V Qt 1 / <f> \ 

M p 2 t (2+pif {m p J ' M p 2 , 2\M P J 

whereas both c| and c T are close to be 1. The tensor ghost can be avoided for ((f>/M p \) p > 0. Under this condition 
the scalar ghost is absent for e 2 = — 1. If £2 = +1, the absence of the scalar ghost requires that 

d1<6. (136) 

For di and d$ of the order of unity we find from Eq. (30) that e2(c?2pids) 2 — 24, where we used the condition (126). 
Hence the dS solution exists only for £2 = +1, in which case d\ ~ 24/(pids) 2 ^ 1- This is incompatible with the 
condition (136). 

These results show that, if the solutions start from the regime |y 2 | <C |x| with e 2 = +1 (i.e. negative kinetic energy), 
the requirement for the avoidance of ghosts at the initial stage is not compatible with the existence of the dS solution 
at late times. 



VII. CONCLUSIONS 

In this paper we have studied the cosmology of generalized Galileon theories based on the Lagrangian (9). For each 
Lagrangian £j (i = 1, • • • , 5) the scalar field <j> is replaced by general scalar functions fi{<j>). The covariant Galileon 
theory satisfies the Galilean symmetry symmetry d^cj) — > d^(j> + in the Minkowski space-time. The extension to 
scalar functions /;(</>) generally breaks this symmetry, but the equations of motion remain at second-order. This is a 
welcome feature to avoid the propagation of the extra ghost degree of freedom. We have also taken into account two 
terms = F(<fi)R and £7 = £(</>)<? that give rise to second-order equations and vanish in the Minkowski space-time. 




(135) 
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In the flat FLRW cosmological background we have derived the equations of motion (11)-(13) for the general 
Lagrangian (9). If we demand the existence of dS solutions, the functions F(<f>), fi (</>), and are restricted to be 
either in the form (17) or (26). The former corresponds to the covariant Galileon theory with constant F, respecting 
the Galilean symmetry in the Minkowski space-time. The latter can be regarded as a kind of scalar-tensor theories 
in which F is field-dependent. 

In the presence of two perfect fluids we have also derived conditions for the avoidance of ghosts and Laplacian 
instabilities associated with scalar and tensor perturbations. The no-ghost conditions (36) and (37) are automatically 
satisfied for the perfect fluids of radiation and non-relativistic matter. Then the no-ghost condition of the scalar mode 
is given by Eq. (38), whereas the ghost is absent for the tensor mode under the condition (42). The stability conditions 
for scalar and tensor perturbations are given, respectively, by Eqs. (44) and (47). We have applied these results to 
two theories having dS solutions. For the theory with constant F the dS solutions are always classically stable against 
homogeneous perturbations, whereas for the theory with non-constant F they are stable under the condition (58). 

We have carried out detailed analysis for the cosmological dynamics of the covariant Galileon theory with constant 
F. Introducing the dimensionless variables r\, and f2 r together with the constants a and /3, it is possible to express 
autonomous equations as well as physical quantities (both background and perturbations) in terms of those variables 
in a convenient form. In particular we showed the existence of an interesting tracker solution r± = 1, along which the 
field velocity evolves as <f> cx 1/H. On this tracker all the non- linear field Lagrangians contribute to the field energy 
density with the similar order, such that any of these terms cannot be neglected. Moreover the cosmological dynamics 
along n = 1 does not depend on the parameters a and j3, see Eqs. (87) and (88). The solutions with different initial 
conditions converge to a common trajectory, depending on the epoch at which they reach the regime n ~ 1. 

Along the tracker solution the dark energy equation of state is given by Eq. (93), which exhibits peculiar evolution: 
u>de = —7/3 (radiation era), wde = —2 (matter era), and wde = —1 (dS era). Since we have derived analytic 
formulas for wde as well as T2 and f2 r in terms of the scale factor a, this will be convenient to confront the Galileon 
theory with supernovae observations. 

Although the background dynamics on the tracker does not depend on the parameters a and /3, the conditions 
for the avoidance of ghosts and Laplacian instabilities do. In Fig. 1 we showed the viable parameter space in the 
{a, /3) plane constrained by the no-ghost and stability conditions along r% = 1. If the solutions start from the regime 
r\ <C 1, we also require the condition /3 > to avoid the scalar ghost. In this case the tensor mode becomes slightly 
super-luminal. In the Minkowski space-time the only solution to the field equation in Galileon theory with d-2 =/= 
corresponds to <p = 0, so that the super-luminal propagation is absent. 

We have also studied the cosmology based on the theories with non-constant F((f>) having de Sitter solutions at late 
times. For the initial conditions with y 2 > i we require that d$px > in the early cosmological epoch. If px > 0, 
there are some viable cosmological trajectories along which the solutions fulfill all the required conditions by today. 
Such an example is given in Figs. 6-9, along which the quantity px remains to be positive. In this case the scalar 
perturbation is subject to the Laplacian instability at the de Sitter fixed point in future (c 2 = — px^/TI < 0). If 
px < initially, we find that the violations of the conditions c| > or Qs > typically occur by today. For the initial 
conditions with y 2 -C x the condition for the avoidance of ghosts in the early cosmological epoch is not compatible 
with the existence of the late-time de Sitter solutions. 

The field-derivative couplings with the Ricci scalar R and the Einstein tensor G up appearing in the terms £4 and £5 
can lead to imprints on the dynamics of matter density perturbations through the change of the effective gravitational 
coupling. It will be of interest to study the evolution of perturbations in detail in order to discriminate between the 
generalized Galileon model and other dark energy models. 
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